.•ur 



- ■ 





NAVAL POSTGRADUATE SCHOOL 
Monterey , California 




THESiS 



6(W2 3 



Investigation of a Heat Driven Thermoacoustic 


Prime Mover 




by 




Lin, Hsiao-Tseng 








December 1989 




Thesis Advisor: 


A. A. Atchley 


Co-advisor 


T. J. Hofler 



Approved for public release; distribution unlimited. 



TP-47274 



Unclassified 

Security' Classification of this page 



REPORT DOCUMENTATION PAGE 



la Report Security Classification Unclassified 


lb Restrictive Markings 


2a Security Classification Authority 


3 Distribution Availability of Report 

Approved for public release; distribution is unlimited. 


2b Declassification/Downgrading Schedule 


4 Performing Organization Report Number(s) 


5 Monitoring Organization Report Number(s) 


6a Name of Performing Organization 

Naval Postgraduate School 


6b Office Symbol 
(If Applicable) 61 


7a Name of Monitoring Organization 

Naval Postgraduate School 


6c Address (city, state, and ZIP code ) 

Monterey, CA 93943-5000 


7 b Address (city, state, and ZIP code) 

Monterey, CA 93943-5000 


8 a Name of Fundmg/Sponsormg Organization 


8b Office Symbol 
(If Applicable) 


9 Procurement Instrument Identification Number 


8 c Address (city, state, and ZJP code) 


1 0 Source of Funding Numbers 


Program Element Number | Project No | Ta»k No | Work Umt Accession No 



1 1 Title Jbiclude Security Classification) INVESTIGATION OF A HEAT DRIVEN THERMOACOUSTIC PRIME 
MOVER 



12 Personal Authors) Lin, Hsiao-Tseng 


13a Type of Report 

Master's Thesis 


13b Time Covered 
From To 


14 Date of Report (year, month,day) 

December 1989 


1 5 Page Count 

75 


16 Supplementary Notation 


Tie views expressed in this thesis are those of the author and do not reflect the official 



policy or position of the Dei 


partment of Defense or the U.S. Government. 


1 7 Cosati Codes 


1 8 Subject Terms (continue on reverse if necessary and identify by block number) 

Acoustics, Thermoacoustics, Thermoacoustic Heat Transport 


Field 


Group 


Subgroup 



1 9 Abstract (continue on reverse if necessary and identify by block number) 

The goal of this thesis is to investigate the work output of a heat driven thermoacoustic prime mover. The 
experimental approach was to measure the frequency response of both a simple resonant tube and a prime mover 
for a variety of values of mean gas pressure and applied temperature difference across the prime mover stack. A 
least squares fit to the frequency response yields the quality factor which can be compared to predictions based on 
a short stack, boundary layer approximation theory by Swift [J. Acoust. Soc. Am. 84, 1 145-1180 (1988)]. The 
results are reported of measurements made on the lowest three modes of the prime mover in helium for mean gas 
pressures between approximately 170 kPa and 500 kPa and the applied temperature differences between zero and 
onset. The signal waveforms of the sound generated by the prime mover above onset at a mean gas pressure of 
308 kPa are also reported. Results for the resonant tube have at most 3% difference with theory. For the prime 
mover, the measurements generally agree with predictions for the fundamental mode except close to onset. This 
agreement between measured and predicted results worsens with decreasing mean gas pressure. Agreement is 
poor for the second and third modes for all pressures used. Finally, the sound generated by the prime mover 
above onset is highly distorted, and the distortion becomes more serve as the temperature difference increases. 
The peak positive pressure amplitude of this signal at temperature difference of 325 °C, 368 °C and 453 °C are 
1.1%, 4.4% and 7.9% of mean gas pressure, respectively. 



20 Distribution/Availability of Abstract 

|~x| unclassified/unlimited same as report PJ DT1C users 


21 Abstract Security ClassificaUon 

Unclassified 


22a Name of Responsible Individual 

Anthony A. Atchley 


22b Telephone (Include Area code) 

(408) 646-2848 


22c Office Symbol 

61 Av 



DD FORM 1473, 84 MAR 83 APR edition may be used until exhausted security classification of this page 



All other editions are obsolete Unclassified 



1 



Approved for public release; distribution is unlimited. 



Investigation of a Heat Driven Thermoacoustic 

Prime Mover 

by 



Lin, Hsiao-Tseng 
Captain, Taiwan Army 

B.S., Chung Cheng Institute of Technology in Taiwan, 1984 



Submitted in partial fulfillment of the requirements 
for the degree of 

MASTER OF SCIENCE IN ENGINEERING ACOUSTICS 

from the 



NAVAL POSTGRADUATE SCHOOL 
December 1989 



ABSTRACT 



The goal of this thesis is to investigate the work output of a heat driven 
thermoacoustic prime mover. The experimental approach was to measure the 
frequency response of both a simple resonant tube and a prime mover for a 
variety of values of mean gas pressure and applied temperature difference across 
the prime mover stack. A least squares fit to the frequency response yields the 
quality factor which can be compared to predictions based on a short stack, 
boundary layer approximation theory by Swift [J. Acoust. Soc. Am. 84, 1145- 
1180 (1988)]. The results are reported of measurements made on the lowest 
three modes of the prime mover in helium for mean gas pressures between 
approximately 170 kPa and 500 kPa and the applied temperature differences 
between zero and onset. The signal waveforms of the sound generated by the 
prime mover above onset at a mean gas pressure of 308 kPa are also reported. 
Results for the resonant tube have at most 3% difference with theory. For the 
prime mover, the measurements generally agree with predictions for the 
fundamental mode except close to onset. This agreement between measured and 
predicted results worsens with decreasing mean gas pressure. Agreement is poor 
for the second and third modes for all pressures used. Finally, the sound 
generated by the prime mover above onset is highly distorted, and the distortion 
becomes more serve as the temperature difference increases. The peak positive 
pressure amplitude of this signal at temperature difference of 325 °C, 368 °C and 
453 °C are 1.1%, 4.4% and 7.9% of mean gas pressure, respectively. 
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I. INTRODUCTION 



There are two basic classes of thermoacoustic engines: heat pumps and prime 
movers [Refs. 1-3]. A thermoacoustic heat pump uses a high amplitude acoustic 
standing wave to transport (or pump) heat along, for example, the boundary of a 
plate situated in the standing wave. This acoustically generated heat flow results 
in a thermal gradient being established across the plate. In other words, acoustic 
energy is converted into stored thermal energy, which in turn can be used in a 
number of practical applications. Thermoacoustic prime movers, on the other 
hand, convert stored thermal energy into useful work in the form of sound. 
Such a device is the subject of this thesis. 

Referring to Figure 1, a typical prime mover configuration consists of a 
stack of plates, called the prime mover stack (or, simply, the stack), which is in 
thermal contact with two heat exchangers. In this thesis, the two heat exchangers 
will be called the hot and ambient heat exchangers, because one end of the prime 
mover stack will be held at elevated temperatures while the other end will be 
held at ambient (room) temperature. However, the important quantity is the 
temperature difference across the prime mover stack, not the absolute 
temperatures of either end. Hence, if (what we called) the hot end were to be 
held at ambient temperature, and (what we called) the ambient end held at 
temperatures below ambient, the process described in this thesis would 
fundamentally be the same. An example of an ambient/cold prime mover is the, 
so called, "Hofler tube", described elsewhere [Refs. 3-4]. As also indicated in 
Figure 1, the prime mover stack and heat exchangers are housed within an 
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acoustic resonator. The prime mover stack/heat exchanger/resonator assembly is 
called a (thermoacoustic) prime mover. 

Thermal energy is stored in the prime mover by imposing a temperature 
difference across the stack, the two quantities being proportional. In order for 
the prime mover to produce net positive work, i.e. to produce audible sound, the 
amount of stored energy converted to sound must exceed the amount of acoustic 
energy dissipated by losses in the prime mover. The dominant loss mechanism 
for the type of gases and frequencies of interest here are thermal and viscous 
losses at the resonator walls and the stack and heat exchanger surfaces. The 
prime mover is said to have reached "onset" when the temperature difference 
across the stack is sufficient for the prime mover to generate and sustain 
detectable levels of sound. The use of the word "detectable" is not meant to 
imply that onset is a subtle question of detection thresholds. Our experience is 
that when onset is reached, the observer (and everyone else in the room) knows 
it. 



The Ambient Heat Exchanger 



The Hot Heat Exchanger 






The Ambient End 



The Hot End 



The Prime Mover Stack 



Figure 1 - A typical prime mover configuration 
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The investigation of a prime mover is conveniently divided into two parts: 
below and above onset. The main goal of this paper is to investigate the 
approach to onset of a thermoacoustic prime mover. However, brief mention 
will be made of the behavior above onset. As will be discussed, there are a 
number of added complications above onset. 

In order to investigate the work output of a prime mover, use will be made 
of the fact that the temperature difference across the stack affects the net 
absorption in the prime mover. As the temperature difference increases from 
zero to its value at onset, the net absorption decreases to zero. The absorption 
coefficient can be determined by measuring the quality factor Q of the prime 
mover, which, in turn, can be related to the acoustic power dissipated in the tube. 
The power can be calculated using equations derived by Swift [Ref. 1], which 
express the acoustic power output of a thermoacoustic engine as a function of the 
applied temperature difference across the stack. 

Swift presents a detailed theoretical development in his review' article. 
How'ever, as he show's, the theory results in relatively simple expressions if the 
short stack and boundary layer approximation are made. The theoretical 
predictions presented in this thesis are based on this approximate theory, in order 
to test the ranges of its validity. 
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II. THEORY 



The measurements performed in the prime mover yield its Q as a function of 
the applied temperature difference. However, the quantity most conveniently 
computed from theory is the net acoustic power dissipated by the prime mover. 
The relationship between these two quantities is derived in this chapter. 

The dependence of the acoustic pressure with distance for a plane wave in a 
closed, rigid, cylindrical tube of length L and cross section area S is given by 
[Ref. 5 Equation (9.26)] 



p(x,t) = P A 



cos[k(L-x)] icot 
cos kL 6 



( 1 ) 



where k = k - ia. The input mechanical impedance presented to a driver located 
at x = 0 is therefore given by 



to 



Z mo=-ip m -j7 S cotkL. 



(2) 



In steady state, the power dissipated by the tube is equal to that delivered by the 
driver, which is given by 



E = 



2 2 

PaS R r 



2Z 



(3) 
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where R m o = Re{Z m0 } and Z m0 = I Z m0 I. Under conditions of resonance for a 
rigid tube, sin kL = 0, Equation 3 reduces to 

E = y— ScxL . ( 4 ) 

Pm C 

The absorption coefficient a can be related to the quality factor Q n of the n th 
mode of the tube by a = C0n/2Q n c. (The subscript n will be dropped unless there 
is potential for confusion). Substituting this relationship into Equation 4 and 
solving for 1/Q gives 



Q 



. 2 
4p c 

— E 
P^cdLS 



(5) 



As indicated in Figure 1, a prime mover is comprised of five parts: the 
hot end, the hot heat exchanger, the prime mover stack, the ambient heat 
exchanger, and the ambient end. Accordingly, the power dissipated by the prime 
mover can be separated into five parts as 



E = Eamb + E amb + E stac j c + E hot + Ehot 
end htex htex end 



( 6 ) 



where the five terms represent the power dissipated by the ambient end, the 
ambient heat exchanger, the prime mover stack, the hot heat exchanger, and the 
hot end, respectively. 
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Using results from Swift’s review article [Ref. 1], expressions can be found 
for the five terms on the right-hand side of Equation 6. We begin with Swift's 
expression for the acoustic power generated by a short thermoacoustic engine 
[Ref. 1 Equation 80] 



W = |n5 v Axco cos 2 (kx) — 



Pm C 



(7) 



(y-0 



0 + Ol (,w?)(l-8 v /y.^fly’ 



- 1 - 






5,t h-5 v /y 0 +5 v /2yJ 



tan (kx) 



where £ s — Pm Cp 8 k / Ps c p §s> T — VT / VTcrit, VT —AT / Ax — (Thot ■ T am b)/Ax, 
and 



VT eri ,= ^0)c CO,kx 



Cp 1+/ /y 0 • 



( 8 ) 



It is assumed in deriving this expression that the stack is short compared to 
the radian wavelength (Ax « X/2n), the boundary layer approximation (v 0 » 
5v, y 0 » S K and / » 8 S ) is valid and that AT « T. It should be pointed out that 
this equation differs from Swift's in that we assume the pressure to vary as 
cos(kx), rather than sin(kx). 

Two modifications, involving conventions, must be made to Equation 7 
before it can be applied to our problem. W in Equation 7 is positive when power 
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is generated. However, the convention in Equation 5 is that E is positive when 
power is dissipated. Hence, E = - W. Finally, the sign convention of AT must be 
reversed relative to Swift’s convention because we will measure x from the 
ambient end, whereas Swift measures x from the hot end. In order to 
accommodate this change, a minus sign will be introduced in front of T. With 
these conventions in mind, we arrive at an expression for the acoustic power 
dissipated by a short thermoacoustic engine 



E = -n6„Axco 
4 



Pa 2 
cos (kx) 



Pm C 



(v- 0 

0+0 



1 + 




r 

5 v /y 0 +5 v /2 yo) 




(9) 




f 1 + / / y ^ 2 

— — (kx) 

5 v /y 0 +S v /2yJ 



Again following Swift, the contribution to Equation 6 from the ambient 
and hot ends of the resonator can be found by realizing that the power dissipated 
per unit area e of the resonator can be found by dividing Equation 9 by the 
surface area of the engine FIAx, setting T = 0, / = 0, and assuming that 5 V « y 0 . 
The justification for letting T = 0 is that the temperature is uniform throughout 
these parts of the prime mover. The justification for setting / = 0 is that the (1 + 
//y 0 ) term accounts for the increase in volume velocity upon entering the stack. 
If there is no stack then this term in unnecessary. Finally, the tube radius R 
assumes the role of y 0 in the absence of the stack. For the frequencies, gases, 
mean pressures, and tube radius of interest, 6v/R in on the order of 10' 2 to 10' 3 . 
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These modifications to Equation 9 yield, for the power dissipated per unit area 
of the resonator 



1 s Pa 2 ., , 

e = — 8^ co cos (kx) 

4 2 
Pm C 



— Q-+— (l + //y 0 ) tan 2 (kx) 

0+0 v 



( 10 ) 



Integration of this expression over the surface area of the ambient and hot ends 
of the resonator yields expressions for E^ and E^ 



£amb — 

end 



P a wLSS k ( y- 1) 


*amb sin(2kx amb ) 


4p m c 2 R d+e s ) 


L + kL 



P a 2 0)LS6 v 


^ amb sin(2kXamb) 


, 6 k (y- l) 


4p m c 2 R 


L kL 


L (l+e s ) 



(ii) 



and 



Ehot 

end 



P a 2 coLS8 k (y_i) 


L-Xhot 


sin(2kx hot ) 


4p m c 2 R 0+e s ) 


L 


kL 


P a 2 coLS5 v 

_u 


L - x hot 


sin(2kx hot 




L 


+ kL 


4 Pm C R 





5 k (y- l) 
L (l+e s ) 



( 12 ) 



The thermophysical properties are to be evaluated for either the ambient 
or hot end, depending on which term is to be computed. 

The contributions to Equation 6 from the ambient and hot heat exchangers is 
found by setting T = 0 in Equation 9 (because the temperature is assumed to be 
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uniform across the heat exchangers) and substituting the result into Equation 5. 
The result is 



• P A con5 K Ax 2 

Eamb/hot= COS (kx) 



ht ex 



4 Pm C 



(y-i) + 8 v 



(i + "y,) 



( 1 + £ s) 5 *(l-5 v /y (> +5 v /2y^ 



tan (kx) 



(13) 



where the geometrical properties of the exchangers and the thermophysical 
properties are to be evaluated for either the ambient or hot heat exchangers, 
depending on which terms are to be computed. 

The stack's contribution to Equation 6 is found by substituting Equation 9 
directly into Equation 5 to give 



P A con6..Ax 



' stack ‘ 



cos 2 (kx) 



4 Pm C 



(14) 



x 



: G- 0 




'' . +1 ) 


4- 


O+'/Yo) 2 


- tan 2 (kx) 

1 , 


0*0 




(l-5 v /y 0 +5 v /2y 2 ) j 


5 «i 


(i - 5 v /y 0 +5 v /2yo) 



Because, as discussed later the initial phases of this experiment consist of 
measuring the Q of an empty, rigidly terminated, cylindrical resonator, we give 
an expression for the Q of such a system. The losses in the empty resonator can 
be found by integrating Equation 10 over the surface area of the entire 
resonator. Substitution of the result into Equation 5 yields 
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1 



( 15 ) 



S K (Y- 1) t 5 v i 25 k (Y- 1) 
Qempty" R (1 + £ ) R + L (1 + e ) 

tube v s 7 v s y 



where is identical to Swift's result [Ref. 1 Equation 91]. 

In order to obtain the Q of the prime mover, it will be driven at frequencies 
in the neighborhood of a resonance with a high impedance source and measure 
the amplitude of the output of a high impedance receiver. Both the source and 
receiver are located at the rigid end of the ambient section of the prime mover, 
i.e. at x = 0. The fitting function is found by considering the acoustic impedance 
seen by the driver. The driver used in these measurements is an electret 
transducer, which acts as a constant displacement source, if the driving voltage is 
constant. The volume velocity of the driver Ud is given by 



where and So are the displacement and surface area of the driver, 
respectively. Vd is, therefore, the volume displaced by the driver. The pressure 
is equal to the product of the acoustic impedance Z and the volume velocity. The 
acoustic impedance is related to the mechanical impedance by Z = Z m /S 2 - 
Therefore, the complex pressure amplitude at x = 0 is given by 



U D -$ D S D co-V D co, 



(16) 



2 




(17) 



where use has been made of Equation 2. 
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The magnitude of the pressure at x = 0 is thus 



p = | p(0) 



P m LV D 



0) 



cot kL 
kL 



(18) 



Using the relation k= co/c, koL= n and Q=co 0 /2ac, the factor kL can be related 
to 0) o and Q by 



f to 

kL = 7t — 

W 




(19) 



By further defining A = p m Vd L co 0 2 / S, Equation 19 can be expressed in terms 
of the three parameters A, Q, and co 0 as 



P = A 



co 



CO / 
\ 0 / 



cot 



71 



CO 



CO 

\ o 



. 

1 



2Q 



n 



/ _co i_' 

“o 2Q 



(20) 



A three parameter, least squares fit of this function to the measured data is used 
to obtain Q and co 0 for the analysis discussed below. 

The question may arise as to the validity of using the assumptions associated 
with Equation 1, in particular, that the losses, represented by the imaginary part 
of the complex propagation constant (k = k - ia), are distributed uniformly 
along the length of the tube. Such is certainly not the case in the prime mover. 
The answer lies in the fact that the quantity of interest is the frequency response 
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of the prime mover, not the variation of the acoustic pressure with distance, or 
even the absolute acoustic pressure in the prime mover. Over a small band of 
frequencies, the pressure amplitude at a particular point should depend only on 
the total attenuation in the tube, not on its spatial distribution. When use is made 
of the relation a = 0)/2Qc, a represents the effective absorption coefficient due 
to uniformly distributed losses that would give the same Q as the one actually 
measured. As long as we limit ourselves to finding the frequency response, the 
formalism associated with Equation 1 is valid. 
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III. EXPERIMENT APPARATUS AND PROCEDURE 



The primary goals are to measure the quality factor of the heat driven prime 
mover below onset as a function of the mean gas pressure and the applied 
temperature difference, and the waveform of the sound generated above onset. 
In order to gain confidence in our ability to measure the quality factor, we 
measured the Q of a simple configuration called the resonant tube, shown in 
Figure 2. This resonant tube was then modified by adding the ambient heat 
exchanger and changing the upper portion to one consisting of the prime mover 
stack and the hot heat exchanger. This configuration (Figure 3) is called the heat 
driven prime mover. A final configuration was built to observe the behavior 
above onset. This configuration is virtually identical to the heat driven prime 
mover, except that the electret driver is replaced by a high intensity pressure 
transducer. This chapter is divided into four sections as follows: (A) The 
resonant tube. (B) The heat driven prime mover. (C) Temperature control 
equipment. (D) Instrumentation and procedure. 

A. THE RESONANT TUBE 
1 . The Tube 

The tube is made from two 3.82 cm ID copper tubes, separated by a 
brass ambient heat exchanger container. The two sections of copper tube are 
11.4 cm and 88.1 cm long. The other end of the two sections are fitted with 
flanges which allow them to be soldered to the ambient heat exchanger container. 
The total length of the resonator, including the heat exchanger container, is 
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Elec tret Driver 

n . 



Gosed Cap 







T 
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* 4 — 


■\T H 



Conner Tube avuhhcui n<^u \ 

Microphone Exchanger Container Copper Tube 



Figure 2 - The resonant tube configuration 







Figure 3 - The Heat Driven Prime Mover 
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1.02m (internal dimensions). A copper cap is soldered to the end of the shorter 
section, forming a closed, rigid, termination. There is also a copper end cap for 
the longer section. However, this cap houses an electret driver and a tiny electret 
microphone. The driver and microphone are flush mounted and sealed in the 
end cap with epoxy. The end cap is epoxied to the section of tube. The 
resonator is connected to a gas handling system, through a valve in the ambient 
heat exchanger container, to allow pressurization with helium. A vacuum pump 
evacuated the resonator before filling with helium. A system vent is also 
provided. In order to sense the mean pressure inside the tube, a dial pressure 
gauge and an OMEGA Model PX304-150AV pressure transducer are connected 
to the fill line. 

2. Driver and Microphone 

It is desired that the impedance of the end cap containing the driver and 
microphone be large compared to the highest acoustic impedance of the standing 
wave. This requirement can be accomplished by using electret designs. An 1.9 
cm electret transducer was constructed for use as the acoustic driver. The input 
signal is provided by a Hewlett Packard Model 8904A Multifunction Synthesizer 
and amplified by a Techron Model 7520 power amplifier. A tiny, 0.594 cm 
diameter, Panasonic electret microphone was used to sense the response. The 
output signal is amplified by a signal amplifier with an open loop gain of 100. 

B. THE HEAT DRIVEN PRIME MOVER 

As explained in the previous chapter, the heat driven prime mover consists 
of hot and ambient heat exchangers and a prime mover stack. These devices are 
described in this section. 
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1 . The Heater Section and Hot Heat Exchanger 

The purpose of the hot heat exchanger is to supply heat to the prime 
mover stack. As shown in Figure 4, the hot end of the prime mover consists of 
a nickel heater section and a heat exchanger. The heater section consists of a 
3.82 cm ID, 5.00 cm long, nickel tube. One end of the tube is capped and drilled 
to accommodate a thermocouple probe, used to sense the hot heat exchanger 
temperature. The other end of the heater section contains the heat exchanger 
consisting of 25, 0.051 cm thick, 0.762 cm long nickel plates. The gap between 
each pair of adjacent plates is 0.102 cm. Between each plates there is a 304 
stainless steel spacers with a length of 0.102 cm and a diameter of 0.031 cm. 

2. The Ambient Heat Exchanger 

In order to impose a temperature gradient across the prime mover 
stack, an ambient heat exchanger is employed to maintain one end of the stack at 
a constant ambient temperature. The construction of this heat exchanger is very 
similar to the hot heat exchanger, except it has a length of 1 .02 cm and contains 
25 copper plates. The ambient heat exchanger actually consists of two such 
stacks separated by a 0.15 cm gap as shown in Figure 5. 

3. The Prime Mover Stack 

This stack is the heart of the heat driven prime mover. A temperature 
gradient will be established across this stack to supply the required heat flux. As 
shown in Figure 6, the prime mover stack consists of a cylindrical stainless steel 
shell containing 35, 0.25 cm thick 304 stainless steel plates 3.42 cm long spaced 
by 0.079 cm . 

A complete listing of the specifications of the heat exchangers and the 
prime mover stack is provided in Appendix C. 
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Figure 4 - The Nickel Heater Section 
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Figure 5 - The Ambient Heat Exchanger 
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Figure 6 - Prime Mover Stack Container 
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C. TEMPERATURE CONTROL EQUIPMENT 

Control of the temperature gradient across the prime mover stack is achieved 
by an OMEGA Model 9151 miniature microprocessor temperature controller, a 
HBA Model 202040 heater and a Neslab Model RET-110 constant temperature 
bath. The heater was mounted to surround the nickle heater section. As shown 
in Figure 7, output of the temperature controller was fed to a voltage divider and 
then connected to the Amplitude Modulation input of the function generator. 
Output from the function generator is then connected to the heater through a 
power amplifier. Water is circulated by the constant temperature bath 
through a water jacket which surrounds the ambient heat exchanger. The 
water pipe was also wrapped around the long section of the prime mover to 
maintain a uniform temperature. A type K thermocouple is soldered to the hot 
heat exchanger to sense the temperature of the hot end. Three type E 
thermocouples were glued to the top, middle, and bottom of the long section of 
the prime mover to sense the temperature along that section. The reference 
temperature for the whole system is found by using a 4-wire method to measure 
the resistance of a thermistor mounted on an isothermal aluminum block. It is 
difficult for an operator to monitor temperature using these thermocouples 
because their outputs give voltages rather than actual temperature. Therefore in 
order to monitor actual temperature directly, two type E thermocouples 
thermometers were glued on to the ambient heat exchanger container and the 
bottom of the prime mover. With this temperature control equipment, the 
maximum deviation of the applied temperature difference is 1 °C. 
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FROM POWER AMPLIFIER 







D. INSTRUMENTATION AND PROCEDURE 

As discussed at the beginning of this chapter, the primary goal is to measure 
the quality factor of the prime mover. The measurement technique employed 
was to drive the prime mover at a frequency near resonance and measure the 
steady state amplitude of the microphone output signal. In this manner, the 
steady state frequency response was obtained. The Q is determined by a least 
squares fit of the ideal response to this data. 

The measurements consisted of three phases. The first phase consisted of 
measuring the Q of the empty resonant tube. The purpose of this phase was to 
give confidence in the measurement technique. The second phase consisted of 
measuring the Q of the prime mover as a function of applied temperature 
difference below onset. The third and final phase consisted of measuring the 
waveform of the acoustic signal produced by the prime mover above onset. The 
instrumentation and procedure was for each of these three phases are discussed in 
this section. 

1. The Empty Resonant Tube 

The experimental setup is shown in Figure 8. The data acquisition is 
performed by an PC AT compatible computer. The computer communicated 
with the SR-530 lock-in amplifier, the HP 3457A multimeter and the HP 8904A 
multifunction synthesizer through a GPIB interface. Through the execution of 
the controlling program, a source signal is supplied by the function generator to 
the electret driver. The output voltage from the microphone was amplified by a 
custom-made signal amplifier. All data signals of interest were fed to a HP 
3457A multimeter. The output of the HP 3457A was used as the data input to the 



21 



computer. A KIKUSUI COS6100A oscilloscope was used to monitor the input 
signal to the electret and the output signal from the signal amplifier. 

Before the data acquisition is started, the valve connected to the tube is 
opened and the system is pumped down and refilled with helium to the desired 
pressure three times in order to purge undesired gases. An ASM Model 110 
turbo helium leak detector is then used to ensure no leakage exists. 

When the automatic data acquisition begins, the program initially takes 
voltages from the PX-304 pressure transducer and converts to actual pressure. 
Next, the approximate resonance frequency and the half power bandwidth are 
entered into the computer, which then determines the start and stop frequencies 
and the frequency increment. After the valve is closed the function generator is 
turned on and the data acquisition is started. The program sets the driving 
frequency and measures temperatures, frequency, X and Y (the real and 
imaginary part of the signal amplitude). It also calculates frequency increment 
and repeats the process. The time required to measure the frequency response is 
approximately 5 minutes. 

As discussed earlier, the resonance behavior of the tube as a function of 
frequency is given by 



where P is the acoustic pressure of the tube, A is a scaling constant, and 
kL=7t[co/co 0 +i/(2Q)]. A least squares fit is performed to obtain fo, Q, and A. 




cot (kL) 
kL 



( 21 ) 
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2. The Prime Mover below Onset 

The setup is shown in Figure 7 (the same portion with Figure 8 is not 
shown here). The measurement was performed as follows: First, the 
measurement performed earlier with the simple resonant tube was repeated. 
After that measurement, the heater was mounted around the nickel heater section 
and then encased in insulating material. The temperature controller and the 
constant temperature bath were turned on and set with a specific temperature. 
After thermal equilibrium is achieved, the data acquisition is started. The 
magnitude of the input electrical power was determined by taking the product of 
the voltage and current. The current was measured through a Tektronix Model 
P6021 current probe. 

Since the quality factor increases with temperature, the HP 8904 
synthesizer was replaced by a HP 3325A function generator to get better 
frequency resolution. This measurement was repeated with increasing 
temperature until sound was generated. 

3. The Prime Mover above Onset 

This configuration (Figure 9) was built to observe the behavior of the 
sound generated above onset. The long section of the prime mover was replaced 
by a different section of identical length which has an ENDEVCO Model 8510B- 
5 piezoresistive pressure transducer screwed through the center of the end cap. 
The transducer is housed within another cap which is bolted to the end cap with 
brass flange. A high impedance leak is provided between the resonator and the 
back volume to eliminate dc pressure difference with little effect on acoustic 
pressure differences. 
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In this phase of the experiment, it was desired to measure the steady 
state electrical power input to the heater. This task was accomplished as follows: 

(1) the HP 3314A function generator was set up with a small constant voltage. 

(2) the controller was set to 900 °C to ensure the heater was turned on all the 
time. After reaching heat equilibrium i.e. the temperature of hot end is constant 
, the data acquisition is started. The output signal from the piezoresistive 
pressure transducer was amplified by an AM 502 differential amplifier and 
connected to a HP 3457 A Multimeter and a HP 3561 Dynamic signal analyzer. A 
Tektronix Model 2445A oscilloscope was used to monitor the output signal and 
the dc voltage supplied to the differential amplifier. 
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Figure 8 - Rxperinient setup of the resonant tube 
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Figure 9 - Measurement setup above onset 





IV. RESULTS AND DISCUSSIONS 



The results of the three phases of the experiment will be presented and 
discussed in this chapter. The results for the empty resonator tube will be 
presented first, followed by those for the prime mover below onset. Finally, the 
results for the prime mover above onset will be presented. 

Investigation of the empty resonator tube and the prime mover below onset 
both involve determining the resonance frequency and the quality factor from 
the measured frequency response. An example of the frequency response data is 
shown in Figure 10, which is a graph of the amplitude of the voltage output from 
the electret microphone as a function of drive frequency. Although this 
particular data set is for the prime mover with a 1 98 °C temperature difference 
across the stack, it is representative of the data obtained under all conditions. 
The solid curve represents the ideal response based on the fit to Equation 21. 
Also indicated in the figure are the estimated errors of the fitted resonance 
frequency and Q values. 

A. THE EMPTY RESONANT TUBE 

The Q's of the first three modes of the empty resonator were measured for 
mean gas pressure ranging from approximately 170 to 500 kPa. The results are 
plotted in Figure 11, which is a graph of Q versus the square root of the mean 
gas pressure. The square root of the mean gas pressure was chosen as the 
abscissa because, as seen from Equation 16, Q « 1 / 5 V P m 1/2 . The solid lines 
represent the theoretical Q as computed from Equation 16. In general the 
agreement is quite good, although there is a tendency to slightly over predict the 
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Q, indicating the presence of unaccounted losses. Two measurements were made 
at each pressure. In all but a few cases the two results fall on top of one another. 
The estimated error of the fit is no larger than the size of the symbols 
representing data points. 

An indication of how closely the empty resonator behaves to an ideal, rigidly 
terminated tube is given in Figure 12. The ordinate of this graph is the ratio of 
the measured (fitted) resonance frequency to the ideal frequency f = nc/2L, 
computed using c = 1008 ms* 1 (the sound speed of helium at 20 °C) and L = 1.02 
± 0.01 m. As seen, all of the values are within 2 to 3 % of 1.00. An error of 
approximately 0.4 % can be attributed to the fact that the actual temperature of 
resonator during the measurements was closer to 18 °C than 20 °C. An 
additional source of error may be due to the fact that the total length of 1.02 m 
was arrived at by adding the individual lengths of the ambient end, the ambient 
heat exchanger container, and the hot end measured separately before the 
resonator was assembly by soldering the three parts together. If the three parts 
did not fit together exactly, the actual length of the resonator may have been 
slightly longer than the unassembled length. 

B. THE PRIME MOVER BELOW ONSET 

The results of the measurements of the Q of the prime mover's fundamental 
mode below onset for five different mean gas pressures are presented in Figures 
13 through 17, which show 1/Q as a function of the temperature difference 
across the stack. The solid line represents the theoretical prediction using 
Equation 6 and 11 through 14, which are based on the boundary layer 
assumption and the short stack approximation. The calculation of the 
contributions from the ambient and hot end is straightforward. The 
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thermophysical properties of the gas are computed using either the ambient or 
hot end temperature, whichever is appropriate. In order to calculate the 
contributions from the heat exchangers, x in Equation 13 is set equal to the 
location of the center of the heat exchanger, and again the thermophysical 
properties of the gas are computed using either the ambient or hot end 
temperature, whichever was appropriate. In order to compute the contribution 
from the stack, x in Equation 14 is set equal to the location of the entrance to the 
stack. The temperature gradient across the stack is accounted for in the 
following manner. Other than the propagation constant, the only temperature 
dependent quantities in Equation 14 are the thermal and viscous penetration 
depths. In order to account for their temperature dependency. Equation 14 is 
converted into a differential by replacing Ax by dx. By assuming a linear 
temperature gradient across the stack, the temperature dependence of the 
penetration depths can be converted into an x dependence. The reciprocal of the 
Q is computed by numerically integrating the differential form of Equation 14, 
with x-dependent penetration depths, along the length of the stack. The 
interpretation of this computation is that, in the short stack limit, the acoustic 
parameters in the stack are determined by their values at the stack entrance. The 
only modification necessary is to adjust the particle velocity for the presence of 
the stack and the viscous boundary layer. 

Before the integration of the penetration depth can be accomplished, it is 
necessary to determine its temperature dependence. This dependence is found by 
performing a least squares fit to viscosity and thermal conductivity data obtained 
from an NBS table [Ref. 6]. The result for both the thermal and viscous 
penetration depths in helium is that 8 « T 0 - 85 . 
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There are several important features evident in Figures 13 to 17. The 
overall agreement is good, through the quality of the agreement diminished with 
decreasing mean gas pressure. In each case, there is a tendency to over predict 
the losses at low temperature difference. Finally, at the lower values of gas 
pressure, the data tend to tail off near onset, a feature not predicted by the 
theory. 

In order to understand the discrepancies discussed above, the validity of the 
approximations must be investigated. For the fundamental mode, the short stack 
approximation is well satisfied because the length of the prime mover stack is 
only about 10 % of the radian wavelength. As for the boundary layer 
assumption, the penetration depth for the average stack temperature is plotted as 
a function of the temperature difference in Figures 17 and 18. If y 0 /5 K ^ 15 the 
boundary layer approximation is valid. For the viscous penetration depth, at 
high mean gas pressure the boundary layer assumption is well satisfied, but this 
assumption starts to fail when the mean gas pressure decreases. The boundary 
layer assumption for the thermal penetration depth is not satisfied very well at 
any pressure. 

In order to determine if there is any connection between the failure of the 
boundary layer approximation and the tailing off of the data at low gas 
pressures, the 1/Q data are replotted in Figures 20 through 24. The solid lines in 
these figures represent a linear least squares fit to the low temperature difference 
data, where the viscous boundary layer approximation is valid. These lines fit 
the 500, 376, and 308 kPa data well over the entire range of temperature 
differences, as might be expected from the fact that the boundary layer 
approximation is valid in these cases at almost all temperature differences below 
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onset. However, the 239 and 170 kPa data begin to depart from the straight lines 
at an approximate temperature difference of 200 and 100 °C, respectively. 
Referring to Figure 18, these are the same points where the approximation 
begins to fail. 

The results for the second and third modes of the prime mover are given in 
Figure 25 through 28 with two different pressures of 308 kPa and 170 kPa. The 
agreement is not very good, probably because the short stack approximation is 
violated at high modes. In the second and third modes the prime mover stack is 
20 and 30 % of the radian wavelength. Although the quantitative agreement is 
poor, the theory does predict the correct dependence of 1/Q at AT = 0 and AT at 
onset. 

C. THE PRIME MOVER ABOVE ONSET 

This final series of data contains the results above onset for the prime mover. 
Figures 29 and 30 show the waveform and the spectrum of the sound generated 
by the prime mover above onset. The mean gas pressure is 307 kPa. The 
temperature difference across the stack is 325 °C which is slightly above onset. 
The signal exhibits slight distortion, particularly in the positive cycle (Figure 
29). Figure 30 shows that the difference in spectrum level between the first few 
modes is larger than 15 dB. Figures 31 and 32 show results for a temperature 
difference of 368 °C. The signal is distorted even further and there is more 
energy stored in high modes. At a temperature difference of 453 °C, the signal 
(Figure 33) is distorted sharply in both positive and negative half-cycles. Figure 
34 illustrates the spectrum for the same signal. The difference in spectrum levels 
between for the first few modes has decreased further to less than 6 dB and more 
energy has been spread to higher modes. 
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The signals shown in these figures have very high amplitude. The peak 
positive pressure in Figures 29, 31, and 33 are approximately 3.3, 13.5, and 24.2 
kPa, respectively. These amplitudes correspond to approximately 1.1, 4.4 and 
7.9% of mean gas pressure. 
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Figure 11 - Resonant tube with helium 
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Figure 12 - Resonant tube with helium 
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Figure 13 - Prime mover with helium at 
500 kPa (the fundamental mode) 
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Figure 14 - Prime mover with helium at 
376 kPa (the fundamental mode) 
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Figure 15 - Prime mover with helium at 
308 kPa (the fundamental mode) 
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Figure 16 - Prime mover with helium at 
238 kPa (the fundamental mode) 
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Figure 17 - Prime mover with helium at 
170 kPa (the fundamental mode) 
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Figure 20 ■ Prime mover with helium at 
500 kPa (the fundamental mode) 
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Figure 21 - Prime mover with helium at 
376 kPa (the fundamental mode) 
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Figure 22 - Prime mover with helium at 
308 kPa (the fundamental mode) 
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Figure 23 - Prime mover with helium at 
238 kPa (the fundamental mode) 
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Figure 24 - Prime mover with helium at 
170 kPa (the fundamental mode) 
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Figure 25 - Prime mover with helium at 
308 kPa (the second mode) 
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Figure 26 - Prime mover with helium at 
308 kPa (the third mode) 
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Figure 27 - Prime mover with helium at 
170 kPa (the second mode) 
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Figure 28 - Prime mover with helium at 
170 kPa (the third mode) 
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Figure 29 - Waveform of the sound generated by the prime mover at a 

temperature difference of 325 °C 



52 



START: 0 Sec STOP: 0 rnSec 

X: 2.110G mSec Y: 30.63 mVolt 



RANGE: -31 rJOV STATUS: PAUSED 

MAG RMS: 20 



< 



> 

D 

■o 



> 
E 
N TD 

^ N 

m o 

OJ • 
r \Tj 

*-i to 
I 



M 

X 

o 

o 

cr 



ic 



> 

o 

T3 

cu 

OJ 

0 
E 

1 







> 



in 
o • 

ITi 0) 



X.jeL 



V v 



<> 

> 



m 




UD > 
OJ CD 



> 

O 2 h 
n -o □ 

\ 



LD 

O 



Figure 30 - Spectrum of the sound generated by the prime mover 
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Figure 31 - Waveform of the sound generated by the prime mover at a 
temperature difference of 368 °C 
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Figure 33 - Waveform of the sound generated by the prime mover at a 

temperature difference of 453 °C 
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Figure 34 - Spectrum of the sound generated by the prime mover at a 



temperature difference of 453 °C 
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V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 



A. SUMMARY 

The purpose of this thesis is to investigate the work output of a heat driven 
thermoacoustic prime mover. The experimental approach was to measure the 
frequency response of both a simple resonant tube and a prime mover for a 
variety of values for mean gas pressure and applied temperature difference 
across the prime mover stack. A least squares fit to a the frequency response 
yield the quality factor which can be compared to predictions based on a short 
stack, boundary' layer approximation theory. 

The results are reported of measurements made on the lowest three modes of 
the prime mover in helium for mean gas pressures between approximately 170 
kPa to 500 kPa and the applied temperatures between zero and onset. The signal 
waveform of the sound generated by the prime mover above onset at a mean gas 
pressure of 307 kPa are also reported. 

The overall results can be summarized as follows. Results of the resonant 
tube have at most 3 % difference with theory. For the prime mover, the 
measurements generally agree with predictions for the fundamental mode except 
close to onset where the boundary layer approximation is not satisfied very' well. 
This agreement between measured and predicted results worsens with decreasing 
mean gas pressure. Agreement is poor for the second and third modes for all 
pressures used, the source of the discrepancy may arise from the violation of the 
short stack assumption at higher modes. Finally, the sound generated by the 
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prime mover above onset has been noticeably distorted. The distortion becomes 
more severe as the temperature difference increases. 

B. CONCLUSIONS 

Several conclusions can be drawn from our results. The first conclusion is 
that the method used to determine the Q works well judging by the close 
agreement between theory and measurement for the simple resonant tube. The 
second conclusion is that short stack, boundary layer theory describes the 
fundamental mode of the prime mover below onset fairly well. There is 
generally good agreement with theory for low temperature differences, although 
the agreement worsens somewhat as onset is approached. The agreement also 
worsen as mean gas pressure decreases, probably due to the break down of the 
viscous boundary layer approximation. A third conclusion is that the short stack 
theory does not adequately describe higher modes, although the theory shows the 
same qualitative features as measurement. The final conclusion is that the sound 
generated by the prime mover above onset exhibits a great deal of nonlinear 
distortion. Reasons for this are not fully understood. 

C. RECOMMENDATIONS 

In order to gain better understanding of the thermoacoustic process and to 
determine the source of discrepancy, particularly for higher modes and the 
region close to onset, several recommendations are proposed as follows: 

• Use a theory not limited by short stack and boundary layer 
approximations. 

• Determine whether the nonlinear distortion seen above onset is primarily 
due to the presence of high amplitude waves in the resonator, independent 
of the presence of the stack, or if the stack is the dominant source of the 
distortion, or if both play equal roles. 
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APPENDIX A. PARTIAL LISTING OF THE PHYSICAL 
PROPERTIES OF HELIUM 



TABLE 1: PHYSICAL PREPERTTES OF HELIUM [Ref. 6] 



p 


T 


DEN 


Cv 


Cp 


C 


vise 


COND 


MPA 


K 


KG/M^ 


J/KG-K 


M/S 


PA-S*E+6 


MW/M-K 


0.170 


293 


0.280 


3123 


5197 


1008 


19.6 


152.4 


0.238 


293 


0.392 


3123 


5197 


1008 


19.6 


152.4 


0.308 


293 


0.504 


3123 


5197 


1009 


19.6 


152.5 


0.376 


293 


0.616 


3123 


5197 


1009 


19.6 


152.5 


0.507 


293 


0.833 


3123 


5197 


1009 


19.6 


152.6 
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APPENDIX B. LISTING OF THE PHYSICAL PROPERTIES 
OF MATERIALS FOR THE PRIME MOVER STACK AND 

HEAT EXCHANGERS 



TABLE 2: PHYSICAL PROPERTIES OF MATERIALS [Ref. 7-8] 

MATERIAL DEn(^) SPECIFIC HEAt(- ; - ~ -— ) CONd(-^-) 

v M 3/ v KG*M y V M*K' 



AISI 304 

STAINLESS STEEL 


8027 


451.9 


16.3 


NICKEL 


8890 


443.8 


89.9 


COPPER 


8900 


384.0 


350.0 
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APPENDIX C. LISTING OF THE SPECIFICATIONS OF THE 
COMPOMENTS AND GEOMETRICAL PARAMETERS OF THE 

PRIME MOVER 



TABLE 3: SPECIFICATIONS OF THE COMPONENTS AND GEOMETRICAL 
PARAMETERS OF THE PRIME MOVER 



•PRIME MOVER STACK 


n= 205.8 cm 


Ax = 3.50 cm / = 0.0125 cm 


y 0 = 0.0395 cm 


•HOT HEAT EXCHANGER 


n = 148.3 cm 


Ax = 0.762 cm l = 0.0255 cm 


y 0 = 0.051 cm 


•AMBIENT HEAT EXCHANGER 


n= 83.7 cm 


Ax = 2.032 cm l = 0.0254 cm 


y 0 = 0.051 cm 


•PARAMETERS OF PRIME MOVER 


x hot = 0.950 m 


X amb = 0.883 m X hot 


= 0.946 m 


x stack = 0.925m 


x a £ b = 0.894 m R= 0.01 91m 


L = 1 .0 m 
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